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The biopolymers actin and microtubules are often in an ongoing assembling/disassembling state 
far from thermal equilibrium. Above a critical density this leads to spatially periodic patterns, as 
shown by a scaling argument and in terms of a phenomenological continuum model, that meets also 
Onsager's statistical theory of the nematic-to-isotropic transition in the absence of reaction kinetics. 
This pattern forming process depends much on nonlinear effects and a common linear stability 
analysis of the isotropic distribution of the filaments is often misleading. The wave number of the 
pattern decreases with the assembling/disassembling rate and there is an uncommon discontinuous 
transition between the nematic and the periodic state. 
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Ongoing polymerization and depolymerization of actin 
and microtubule filaments are prominent examples for 
dissipative non-equilibrium phenomena in living cells , 
which are important for many different purposes, such 
as the cell motility and division or morphogenesis. Both 
substances show also an inherent propensity to pattern 
formation and active phenomena 0, IE IE IE Hi - Like 
the famous example for rod-like particles, the Tobacco 
Mosaic Virus (TMV) , also actin and microtubule fila- 
ments may undergo with increasing density a well known 
transition to an orientational order [l^ |lj |l3 , the so- 
called nematic order 0, Q| . By Onsager's seminal work 
^EQjEj .IO] this transition has been traced back to ex- 
cluded volume interactions between the filaments. This 
statistical theory is valid for long filaments of fixed shape 
and infinite lifetime r, and it predicts near the transition 
also a phase separation into domains of isotropically ori- 
ented rods at low density and nematic domains of higher 
rod-density, which has also been observed for actin with 
an almost vanishing kinetics ^E|- For a finite lifetime 
of actin and microtubule filaments, Onsager's statistical 
theory for the nematic order does not apply. Moreover, 
a finite r limits the diffusive transport distance and the 
coarsening during the phase separation close to the orien- 
tational transition to a length scale of about Ijj = ^yOpT, 
with the filament diffusion constant Dp. According to 
our estimate we expect kinetically induced periodic pat- 
terns with a wavelength in the order of 10 /xm. This is 
supported by current experiments . 

The effect of reaction kinetics on a phase separation 
has been investigated for a chemical and a biophysical 
example in Refs. [rRll8| . In both cases the transition to 
periodic patterns is supercritical and its onset as well as 
wavelength follows already from a linear stability analy- 
sis of the respective homogeneous basic state. Instead of 
these two competing states, near the orientational transi- 
tion of filaments one has three different competing states, 
the spatially homogeneous isotropic state, the spatially 
homogeneous nematic one and the spatially inhomoge- 
neous alternation between the isotropic and the nematic 



order. Here, the growth rate of perturbations of the ba- 
sic state takes its maximum also at a finite wavenum- 
ber [1^ l2nl |. but this is not sufficient for a prediction 
of spatially inhomogeneous nonlinear states above the 
orientational transition. Instead of an inhomogeneous 
state, as suggested by the linear perturbation analysis, 
in the nonlinear regime one has an exchange of stabil- 
ity and the spatially homogeneous nematic state is often 
preferred, as described in this work. Therefore the bifur- 
cation picture has to be explored by a nonlinear analysis, 
whereby an uncommon coarsening behavior at the non- 
linear nematic-to-periodic transition has been found. 

The generic scenario near the ordering transition is de- 
scribed in terms of a phenomenological model, which is 
introduced and analyzed at first without the reaction ki- 
netics of the filaments. It is extended in the second part 
by the essential reaction steps as motivated by actin and 
microtubule polymerization. 

Model without reaction kinetics. - In lyotropic liquid 
crystals the nematic order is forced beyond a critical rod 
density pc by excluded volume interactions and the 
resulting local mean orientation of rod-like particles is 
described by the so-called director n(r) (with n = — n) 
|l^. Assuming rods of a single length and a uniform n, 
then the largest eigenvalue A of the nematic order param- 
eter tensor is sufficient for a description of the strength 
of the orientational order. A varies in the range [5,1], 
with A = I in the isotropic state and A = 1 for a uni- 



form rod orientation |14L |21| . It is convenient to use the 
difference A := A — and the dimensionless rod density 
p oc pVe, with Ve — 2bL^ the excluded volume for rods 
of length L and diameter b. For a constant orientation 
n, as we assume in this work, spatial variations of p and 
A also include spatially alternating isotropic and nematic 
ranges. Since for long rods the preferred director orien- 
tation is parallel to the isotropic-nematic interface [2^ . 
we assume spatial variations in the direction perpendic- 
ular to n, which we call the a;-direction. For this context 
we choose the phenomenological model for the conserved 
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density p{x,t) and for the unconserved field A(x,t) 



dtp = D,dl[-\p-Spdlp 



(la) 



{l-p)X--p\^ 
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+Dxdl [(1 - p)A - SxdlX + oaA^ 



(lb) 



In Eq. Ijlb|l . the part without spatial derivatives follows 
from the Smoluchowski equation for rigid rods by a mo- 
ment approximation |2lL |23j and the factor p in front 
of the nonlinear terms reflects the excluded volume in- 
teraction. It determines also the homogeneous solutions 



1 1 I 8 

p^Po^ const., Ao = 0, A± = -±-W9 . (2) 

6 6V Po 

Ao = corresponds to the isotropic rod distribution that 
becomes linearly unstable with respect to ncmatic fluctu- 
ations beyond the critical density p > pc = 1, where they 
grow up to the homogeneously stable upper branch A-f of 
the spatially uniform nematic order in Fig.^). Since the 
isotropic-nematic transition is of first order, both states 
coexist in a range | < p < 1. 

Eq. (|la|l is of the Cahn-Hilliard type Expressing 
its right hand side by a divergence of the current density 
jp{x) = -Dpdxp{x) with p{x) = -Xp - Spd^p + app^, 
it takes the form of a conservation law for the rod-like 
particles. The first (nonlinear) term in Eq. Hla(l . i.e. 
—Dpd^{Xp), destabilizes the spatially homogeneous par- 
ticle density for any finite value of A (A is always posi- 
tive) and mimics therefore Onsagers prediction [ij, 
that the free energy can be reduced, by separating the 
system into ranges of low rod density pi (isotropic) and 
high density pa (nematic). The second term describes 
an isotropic-nematic interface energy and the third term 
limits the modulation amplitudes of the density (see also 
below) . 

Theories involving the distribution function predict be- 
yond the critical pc = 1 an instability of the isotropic dis- 
tribution against inhomogeneous order parameter fluctu- 
ations 0, . In Eq. (|lb|l this is taken into account by 
i9^((l — p)A) and the last two terms limit the wavenum- 
ber and the amplitude of the nonlinear modulations of 
A. For intermediate values of p, Eqs. have stationary 
kink solutions as shown in Fig.^). The densities pi and 
Pa, in the isotropic and nematic range respectively, are 
determined by the two coefficients ap and a\ , as described 
in the following. For stationary kinks as in Fig. ^) the 
particle transport vanishes jp{x) = and p{x) — pi — pa 
is constant. Sufficiently far away from the kink pa and 
Pi are constant too and one obtains the equation 



OpP? = -XaPa + appl , 



(3) 



where Aq = X{pa) = (1 + - 8/ Pa )/6. The rota- 
tional term in Eq. Ijlb|l vanishes in the nematic range 




FIG. 1: Part b) shows the stable (solid) and the unstable 
branch (dashed) of the nematic order parameter X± as a func- 
tion of the rod density p, cf. Eq. J^J. Part c) shows stable 
kink solutions of Eqs. Q interpolating between the nematic 
{p ~ Pa) and the isotropic range (p = pi). The nematic vol- 
ume fraction Vn as a function of p is given in a). Parameters 
of the model: 
2.0,3, = Sx = 



0.1. 



0.1, Do 



Dx 



0.3, a„ 



0.25, 



ax 



For this set, the whole system is in the 



homogeneous nematic state for po > 1.367 (i.e V„ — 1) 

as well as trivially in the isotropic range. Defining 
jx{x) = —Dxdxv{x), since in the isotropic range A^ is 
zero, Vi is zero as well. To prevent a current through the 
interface, the total current, i.e. just jx in the nematic 
region, has to vanish and it follows Va = ~ Q leading 
to 



(1 - pa)Xa + axXl = 



(4) 



As Xa is known, from this equation the anisotropic den- 
sity Pa follows as a function of ax (or vice versa) and pi 
is given with pa via Eq. © as a function of ap (or vice 
versa). Therefore, the two densities pi and Pa may also 
be considered as input parameters that are obtained from 
different approaches as for instance from Refs. |l3tl25i | or 
possibly from experiments. Since pi and pa do not de- 
pend on the system size L, the kink position changes 
with the mean density po as shown in Fig. ^ where the 
nematic volume fraction Vn is given in terms of po . 

Linear stability analysis. - For p > pc & linear stability 
analysis of the homogeneous isotropic state, cf. A = 
and p ~ Po, with respect to small periodic perturbations 
Ai, pi oc exp{at±iqx) gives the wave number dependence 
of the growth rate cr(g) as shown for a set of parameters in 
Fig. lit). This shape of a{q) with a positive value at g = 
and a maximum at a finite value of q is typical for the 
unstable isotropic state and is i n ag reement with similar 
results for microscopic models [l9j. A linear stability 
analysis of the homogeneous nematic state X± in terms 
of microscopic models is rather involved (23|. For our 
phenomenological model, however, the determination of 
cr(g) is a straightforward task and its typical shape at 
the unstable nematic branch is shown in Fig.|2b)- Along 
the dashed part of the curve in Fig. ^p) the homogeneous 
nematic state is linear unstable. Between our result and 
the linear stability described in there is a major 
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difference. In both cases a{q) for the isotropic state takes 
its maximum at a finite value of q and has positive values 
for any p > p^. This is somewhat in contradiction to 
Onsager's statistical theory, where inhomogeneous states 
(via phase separation) are only energetically preferred 
for a rod density below a maximum value p < pa- In 
our model nonlinear effects stabilize the uniform nematic 
state for p > Pa and it is unstable only along the dashed 
line in Fig. ^fl). Simulations of Eqs. confirm that 
inhomogeneous solutions X{x) and p{x) only occur for a 
mean density po smaller than pa- 



a) 




0.2 0.4 0.6 0.8 1.0 q 0.2 0.4 0.6 0.8 Q 0.2 0.4 0.6 0.8 Q 



FIG. 2: In part a) a{q) is shown for periodic perturbations of 
the isotropic state with po > pc and in part b) of the nematic 
state without kinetics and at the unstable branch in Fig. 0d) . 
In part c) (j{q) is given for perturbations of the nematic state 
in the case with reaction kinetics and at the unstable branch 
in Fig. In part a) and b) the parameters are po = 1.05 
and E = and in c) E = 0.003 and s = O.OIE. The other 
parameters are as in Fig. 

Reaction kinetics drives pattern. - In cells and in vitro 
actin and microtubule filaments are usually out of equi- 
librium and, due to an ongoing assembly/disassembly 
reaction, filaments have a finite lifetime t. This re- 
action kinetics leads to a stationary length distribu- 
tion of the filaments or even to oscillatory poly- 
merization 0, ■ During the phase separation at 
the isotropic-to-nematic transition, filaments are trans- 
ported, but only over a lifetime-dependent distance of 
about Id = ^DpT. Since the lifetime of filaments is a 
constant, much more subunits are released in the nematic 
range with a high density pa than in the isotropic range 
with a low density pi. However, due to a much larger 
diffusion constant, the subunits are redistributed quickly, 
leading to a nearly homogeneous subunit density m{x). 
Thus the number of nucleated filaments per unit time, 
which depends on m(x\ is weakly varying too. By this 
qualitative reasoning one expects a steady net transport 
of subunits from the nematic to the isotropic range and in 
the opposite direction a transport of filaments, whereby 
the latter one is limited to distances of the order of I u or 
smaller. This length restriction causes, instead of a large 
scale phase separation, a spatially periodic pattern with 
a wavelength in the order of Zd. 

Along this qualitative reasoning the complexity of the 
biochemical reaction steps, involved during the assem- 
bly/disassembly of actin or microtubules, is not crucial 
for this wavelength limitation. For instance, actin and 
microtubules are usually met with a polydisperse length 
distribution. Since the slowest kinetic step and the small 
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FIG. 3; Part a) and b) show periodic solutions p(x) and \{x) 
of Eqs. Q for po ~ 0.8 (dotted), po = 0.95 (dashed) and 
Po = 1.056 (sohd), respectively. Part d) displays the stable 
homogeneous nematic branch (solid) and the unstable ones 
(dashed). In part c) the existence range of the stable (solid) 
and unstable (dashed) periodic patterns with the modulation 
amplitude A of A(a;) are given. Parameters are as in Fig. 
with L>m = 10, 7 = 100, E = 0.003, s = O.OIS. 

diffusion constant of the long filaments will govern the 
limitation, we discard the polydispersity and assume for 
the sake of simplicity that all filaments are of the same 
length. With a decay S — and a nucleation rate s of 
the filaments and a diffusion constant Dm of the subunits 
one ends up finally with the three equations 

dtp = Dpdl [-Ap - 5pdlp + app^] + sm- Sp, (5a) 



dt\ = -Dr 



(l-p)A-^pA^ + ^pA^ 



+D),dl [(1 - p)A - SxdlX + oaA^] " SA , (5b) 
7sm + 7SP . (5c) 



The constant 7 is a measure for the number of subunits 
per rod-like particle. In the nematic range oriented fil- 
aments are lost, but new ones are nucleated everywhere 
with an arbitrary orientation which have to relax to the 
local mean orientation by rotational diffusion. Accord- 
ingly, there is only a decay term in Eq. (|5b|) which can be 
justified also microscopically |2^. This reaction kineti- 
cally caused partial loss of the orientational order leads 
to a higher critical density for the isotropic-to-nematic 
transition pc = 1 -l- H/ D^- In the presence of reactive 
steps the rates s and S determine the mean rod den- 
sity in terms of the monomer density p — sYi^^m. The 
spatially homogeneous solutions of Eqs.Q are 



P = Po, Ao=0, A. = i4^9-A(l-K|),(6) 

where A± (p) are depicted as a function of the mean den- 
sity p in Fig.Oi). In a certain parameter range A+ may 
become unstable with respect to periodic perturbations 
Ai, pi oc exp(f7t ± iqx) and the wave number depen- 
dence of the growth rate u{q) has a typical shape as 
in Fig. 13;). Compared to the case without kinetics as 
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in Fig. I^t), long wavelength perturbations are now sup- 
pressed and only perturbations with a finite wave number 
grow. The parameter range of the unstable homogeneous 
A_|_ branch is indicated by the dashed and the stable one 
by the solid line in Fig. Ol). The unstable range of the 
nematic branch decreases with increasing coefficients 6p 
and 6x and with increasing and moderate values of the 
decay rate E. Periodic states occurring at the unstable 
nematic branch are shown in Fig.OJi) and b) for three dif- 
ferent mean densities pq. The maxima of p{x) come close 
to Pa, cf. Fig.^ but the minima of are considerably 
larger than pi in the case without kinetics. The solid line 
in Fig. IHt) indicates the range where the nonlinear peri- 
odic state is in coexistence with the homogeneous states. 
For p > p*, along the dashed line in Fig. the peri- 
odic pattern becomes increasingly anharmonic, plateaus 
P ^ Pa spread out and the valleys of low filament density 
in between become less and narrower by approaching the 
upper end of the dashed curve in Fig.^t), a behavior that 
is rather uncommon |2^. At the left end of this curve 
the state remains periodic, the wavelength increases and 
the valleys spread out. 

Conclusions.- A reaction kinetically driven pattern 
forming process is predicted near the isotropic-nematic 
(I-N) transition, which is supported by recent experi- 
ments on actin polymerization (Si]- A phenomenological 
continuum model is introduced that reproduces the first 
order I-N transition in lyotropic liquid crystals, includ- 
ing the phase separation in its neighborhood and being in 
agreement with the statistical theory of Onsagcr. Peri- 
odic solutions arise due to a finite lifetime r and a nucle- 
ation rate of the filaments. However, the correlation be- 
tween a finite t and the occurrence of periodic patterns is 
independent of the details of the model. Beyond a critical 
density pc the isotropic orientation of the filaments and 
below a certain p* the uniform nematic state becomes un- 
stable against inhomogeneous perturbations. Hence for 
p < p* periodic states are favored and for our model the 
wavenumber varies as q cx r^*^'^^. However, pc increases 
and p* decreases with E = r^^ and it may happen that 
p* < Pc holds, i.e. in a certain parameter range kinetics 
favors the uniform nematic state. Therefore, inhomoge- 
neous states as predicted by a linear stability analysis 
of the isotropic state [l^ |20| may be irrelevant due to 
nonlinear effects. This pattern formation process near 
the I-N transition is expected to be generic and may also 
apply to situations with different transport and filament 
accumulation mechanisms such as in the case of filament 
bundling [29I . Isol Isif . Even though our description is 
very simplified and restricted to one spatial dimension, 
we expect that the basic physical mechanism also plays 
a crucial role for situations with polydisperse filament 
distributions 0, El..22J, including living cells. Polydis- 
persity favors periodic patterns and together with higher 
spatial dimensions this will give rise to an even larger 
variety of phenomena, to which forthcoming works are 



devoted. 
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